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The  mechanisms  of  formation  and  destruction  of  NO  in  MILD  combustion  of  CHJH2  fuels 
blends  are  investigated  both  experimentally  and  numerically.  Experiments  are  carried  out 
at  a  lab-scale  furnace  with  the  mass  fraction  of  hydrogen  in  fuel  ranging  from  0%  to  15%; 
furnace  temperature,  extracted  heat  and  exhaust  NOx  emissions  are  measured.  Detailed 
chemical  kinetics  calculations  utilizing  computational  fluid  dynamics  (CFD)  and  well- 
stirred  reactor  (WSR)  are  performed  to  better  analyze  and  isolate  the  different 
mechanisms. 

When  the  MILD  combustion  of  the  CH4/H2  fuel  is  established  in  experiments,  the 
thermal  field  is  quasi  uniform  and  the  high  temperature  zone  is  located  at  the  junction  of 
the  fuel  and  air  jets.  As  the  mass  fraction  of  hydrogen  in  fuel  is  increased  from  5.7%  to 
14.4%,  although  the  furnace  average  temperature  is  increased,  the  NOx  emission  remains 
unchanged.  This  cannot  be  explained  by  the  thermal  NO  mechanism.  CFD  and  WSR  sim¬ 
ulations  both  suggest  that,  when  equivalence  ratio  <0.8,  the  N20-intermediate  route 
controls  the  NO  formation  and  the  NO-reburning  reaction  is  also  strong.  With  the  hydrogen 
addition,  the  importance  of  the  NNH  route  is  increased  but  that  of  the  prompt  route  is 
decreased,  consequently  non-affecting  the  NOx  emission  as  measured. 

Chemical  kinetics  calculations  indicate  that  the  conversion  from  NO  to  N02  becomes 
significant  and  thus  the  relative  importance  of  N02  is  increased  in  the  total  NOx  emission 
under  low  temperature  MILD  conditions.  As  the  reactor  temperature  is  increased  from 
1100  K  to  1600  K,  the  importance  of  N20  route  decreases  while  that  of  thermal  route  in¬ 
creases.  In  contrast,  as  the  initial  mass  fraction  of  oxygen  is  increased  from  3%  to  9%,  the 
importance  of  N20  route  increases  but  that  of  the  prompt  and  NNH  routes  decreases. 
Likewise,  as  the  equivalence  ratio  increases,  the  NO-reburning  reaction  becomes  strong. 
Worth  noting  is  that  the  N20-intermediate  route  controls  the  NO  production  under  fuel 
lean  conditions  whereas  the  prompt  route  is  dominant  in  rich  ones. 

Copyright  ©  2014,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights 

reserved. 


*  Corresponding  author.  State  Key  Laboratory  of  Turbulence  and  Complex  Systems,  Department  of  Energy  &  Resources  Engineering, 
College  of  Engineering,  Peking  University,  Beijing,  100871,  China.  Tel./fax:  +86  01062767074. 

E-mail  address:  jcmi@coe.pku.edu.cn  (J.  Mi). 
http://dx.doi.Org/10.1016/j.ijhydene.2014.09.050 

0360-3199/Copyright  ©  2014,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights  reserved. 


19188 


INTERNATIONAL  JOURNAL  OF  HYDROGEN  ENERGY  39  (2OI4)  I9187-I9203 


Nomenclature 

Symbols 

A i  pre-exponential  factor 

Q  the  volume  fraction  constant  of  eddy  dissipation 

concept  model  (=2.1377) 

Cr  the  time  scale  constant  of  eddy  dissipation 

concept  model 

Da  diameter  of  the  air  nozzle  exit,  mm 

Df  diameter  of  the  fuel  nozzle  exit,  mm 

Kv  gas  recirculation  rate 

p  pressure,  atm 

Qextracted  extracted  heat,  kW 
Qinput  thermal  input,  kW 

T  temperature,  K 

f  average  average  temperature  of  the  furnace  centerline,  K 

"f exhaust  exhaust  temperature,  K 
'f mean  mean  temperature  of  the  entire  furnace,  K 
Twsr  temperature  of  the  well-stirred  reactor,  K 


XNox  NOx  emission  in  volume,  ppm 
YCh4  initial  mass  fraction  of  CH4  in  the  well-stirred 
reactor,  % 

YCo2  mass  fraction  of  C02,  % 

YfH  initial  mass  fraction  of  H2  in  fuel,  % 

Yh2  initial  mass  fraction  of  H2  in  the  well-stirred 

reactor,  % 

Yj  the  mass  fraction  for  the  jth  species,  % 

YN2  initial  mass  fraction  of  N2  in  the  well-stirred 
reactor,  % 

Yo2  initial  mass  fraction  of  02  in  the  well-stirred 
reactor,  % 

Yg2  initial  mass  fraction  of  02  in  oxidant,  % 

Greek  letters 

£  characteristic  length  fraction  of  fine  scales,  m 

r  residence  time,  s 

0  equivalence  ratio 


Introduction 

Binary  fuel  blends  are  widely  used  in  developing  high  effi¬ 
ciency  and  low  emission  combustion  devices  [1],  Adding 
hydrogen  to  methane  (or  natural  gas)  can  promote  ignition 
and  enhance  flame  stability  [2—5],  and  offer  definite  advan¬ 
tages  in  applications  such  as  combustion  furnaces  and  en¬ 
gines  [6,7].  However,  the  hydrogen  addition  generally 
increases  the  flame  temperature  and  thus  yields  more  NOx 
emissions.  To  solve  this  problem,  one  of  the  effective  methods 
is  to  operate  the  combustion  in  moderate  or  intense  low  ox¬ 
ygen  dilution  (MILD)  condition  to  suppress  NOx  emissions 
[8-10]. 

The  MILD  combustion  is  achieved  by  strong  exhaust  gases 
recirculation,  where  reactants  are  intensely  diluted  and  hence 
reactions  occur  volumetrically  without  visible  flame-front. 
Relative  to  the  conventional  counterpart,  the  MILD  combus¬ 
tion  increases  the  in-furnace  thermal  uniformity  and  effi¬ 
ciency,  and  simultaneously  suppresses  NOx  emission  [11—14]. 
Although  considerable  work  has  been  published  on  MILD 
combustion  [14—28],  most  of  the  previous  investigations  have 
focused  on  firing  hydrocarbon  gaseous  fuels  or  pulverized 
coal.  To  understand  the  characteristics  of  the  MILD  combus¬ 
tion  of  hydrogen  containing  fuels,  both  the  experimental  and 
numerical  investigations  are  essential. 


In  fable  1  a  summary  of  published  experimental  work  on 
the  MILD  combustion  of  the  CHyH2  fuel  in  furnaces  is  pre¬ 
sented.  Derudi  et  al.  [29,30]  showed  that,  relative  to  burning 
pure  CH4,  the  MILD  combustion  of  the  CH4/H2  fuel  requires 
higher  jet  velocity  to  establish,  but  it  can  operate  at  lower 
average  furnace  temperatures.  They  also  found  that  the 
addition  of  hydrogen  leads  to  complete  oxidation  of  the  hy¬ 
drocarbon  under  MILD  conditions.  Parente  et  al.  [31]  carried 
out  numerical  and  experimental  investigations  of  the  MILD 
combustion  of  the  CH4/H2  fuel  (with  hydrogen  content  up  to 
20%  by  wt.).  For  their  MILD  combustion  system,  they  discov¬ 
ered  that  the  influence  of  molecular  diffusion  should  be 
considered  in  the  numerical  prediction  of  H2  distribution,  but 
its  effects  on  the  temperature  field  and  major  species  are 
negligible.  They  also  developed  a  simple  NO  formation 
mechanism  based  on  the  thermal  and  prompt  routes  to  pre¬ 
dict  NO  emissions.  The  MILD  combustion  in  their  experiments 
notably  operated  in  relatively  high  temperatures  (>1800  K). 
Parente  et  al.  [32]  further  developed  a  simplified  approach  to 
predict  NO  formation  in  the  MILD  combustion  of  CH4/H2 
mixtures.  They  found  that  the  NNH  and  N20-intermediate 
routes  are  crucial  when  the  MILD  combustion  occurs  at  low 
temperatures  and  hydrogen  is  added  to  the  fuel.  Arghode  et  al. 
[33]  found  that  the  addition  of  H2  results  in  a  slight  increase  of 
the  NO  emission  for  the  same  equivalence  ratio  possibly  due 
to  a  higher  flame  temperature.  Galletti  et  al.  [34]  pointed  out 


Table  1  -  Summary  of  experimental  investigations  on  the  MILD  combustion  of  the  CH4/H2  fuel  in  furnaces. 

Fuel 

Co* 

CM 

£ 

CN 

D(K) 

P  (kW) 

NOx  emissions  (ppm) 

CO  emissions  (ppm) 

Reference 

CH4/H2  mixture 

5,  15 

>1173 

0.2— 0.3 

<30 

<50 

[29] 

CH4/H2  mixture 

15 

>1173 

CO 

d 

1 

CM 

d 

<30 

<50 

[30] 

CH4/H2  mixture 

5.5 

800-900 

13 

59-105 

N/A 

[31] 

CH4/H2  mixture 

50 

N/A 

10 

25-45 

N/A 

[32] 

CH4/H2  mixture 

8 

300 

6.25 

0-30 

20-600 

[33] 

CH4/H2  mixture 

15 

1223,  1312 

CO 

d 

1 

CM 

d 

17 

<50 

[34] 

CH4/H2  mixture 

0-100 

298,  858 

20 

8-16,  26-42 

3-60,  <10 

[35] 
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that,  for  the  MILD  combustion  operating  in  1300—1400  K,  the 
NNH  and  N20-intermediate  routes  are  the  dominant  forma¬ 
tion  pathways  for  the  MILD  combustion  of  CH4/H2  mixtures. 
Ayoub  et  al.  [35]  carried  out  MILD  combustion  experiments  in 
a  laboratory-scale  furnace  (20  kW)  and  found  that  MILD 
combustion  can  be  established  for  CH4/H2  mixtures  whatever 
the  mass  fraction  of  H2  and  even  when  pure  H2  was  used.  In 
their  experiment,  with  the  addition  of  H2,  the  NOx  emission 
was  decreased  due  to  the  decrease  of  prompt  NO  route.  Yu 
et  al.  [36]  performed  a  well-stirred  reactor  (WSR)  network 
simulation  on  the  MILD  combustion  of  CH4/H2  mixtures.  They 
found  that  CO  decreases  from  48  ppm  to  zero  when  the  vol¬ 
ume  fraction  of  H2  is  increased  from  40%  to  100%,  but  the  NOx 
emission  is  not  changed  significantly.  They  also  discovered 
that  the  NOx  and  CO  emissions  do  not  depend  on  the  exhaust 
gas  recirculation  in  the  highly  diluted  conditions. 

In  addition  to  studies  in  furnaces,  the  investigation  of  MILD 
combustion  was  also  performed  using  jet-in-hot-coflow 
burner  which  emulates  the  MILD  combustion  regime  locally 
in  furnace  [23,37,38].  Medwell  et  al.  [39]  carried  out  laser 
diagnostic  experiments  using  a  jet-in-hot-coflow  burner 
operating  at  MILD  conditions.  Their  experiments  demon¬ 
strated  that  H2  addition  to  the  primary  fuel  can  stabilize  the 
flame,  but  the  structure  of  the  reaction  zone  is  not  varied 
significantly.  Mardani  et  al.  [40]  numerically  investigated  the 
CH4/H2  jet  in  hot  and  diluted  coflow  air  flame  under  the  MILD 
condition.  These  authors  revealed  that  the  hydrogen  addition 
to  methane  increases  reaction  rates.  Mardani  et  al.  [41]  also 
pointed  out  that  the  NNH  and  N20-intermediate  routes  are 
most  important  pathways  in  NO  formation  under  MILD  con¬ 
ditions.  Moreover,  they  discovered  that,  as  the  hydrogen  level 
is  decreased,  the  importance  of  the  NNH  and  N20-interme- 
diate  routes  decreases.  The  investigation  of  Gao  et  al.  [42] 
indicated  that  the  NNH  and  prompt  routes  play  a  significant 
role  in  the  NO  formation  under  the  MILD  conditions  of  CH4/H2 
mixtures  using  the  jet-in-hot-coflow  burner.  It  is  worth  noting 
that  the  mass  fraction  of  the  oxygen  in  the  hot  coflow  is  not 
low  (9%).  Their  results  also  revealed  that  the  NO  formation  is 
suppressed  as  the  mass  fraction  of  hydrogen  is  decreased. 
Sepman  et  al.  [43]  investigated  the  effect  of  hydrogen  addition 
on  the  flame  structure  of  MILD  condition  using  a  laminar-jet- 
in-hot-coflow  burner.  They  found  that  the  addition  of 
hydrogen  substantially  decreases  the  flame  height  (-25%)  but 
only  slightly  affects  the  maximal  flame  temperature  and  the 
thickness  of  combustion  zone.  Afarin  et  al.  [44]  examined  the 
effect  of  hydrogen  addition  on  the  CH4/H2  flame  structure 
under  the  MILD  condition  using  large  eddy  simulation.  These 
investigators  concluded  that  the  addition  of  hydrogen  to  fuel 
not  only  increases  both  the  peak  temperature  and  the  flame 
thickness  but  also  reduces  the  fluctuation  of  the  OH  radical 
and  the  diffusion  of  the  un-burnt  fuel  in  the  flame  front. 

The  previous  studies  cited  above  have  preliminarily 
investigated  the  NO  formation  mechanisms  under  the  MILD 
combustion  of  CH4/H2  fuel  blends.  However,  no  systematical 
work  has  been  conducted  that  examines  the  effects  of  various 
operational  parameters  (i.e.,  the  hydrogen  addition,  temper¬ 
ature,  oxygen  level,  and  equivalence  ratio)  on  the  formation 
routes  of  NO.  The  present  study  is  our  effort  to  help  address 
this  issue.  Two  specific  objectives  are  designated:  (1)  to 
experimentally  investigate  the  effect  of  hydrogen  addition  on 


the  performance  of  MILD  combustion,  and  (2)  to  numerically 
examine  the  NO  mechanism  of  the  MILD  combustion  of  CWJ 
H2  mixtures  systematically.  For  the  first  objective,  experi¬ 
ments  are  performed  at  the  mass  fraction  of  hydrogen  in  fuel 
ranging  from  0%  to  15%.  The  flow  field,  furnace  temperature, 
extracted  heat  and  exhaust  NOx  emissions  are  obtained.  For 
the  second  objective,  the  NO  mechanism  is  investigated  by 
CFD  and  WSR  simulations.  In  particular,  effects  of  hydrogen 
addition,  reactor  temperature,  oxygen  level,  and  equivalence 
ratio  on  the  NO  formation  and  reduction  mechanisms  under 
the  MILD  condition  of  CH4/H2  mixtures  are  investigated. 
Finally,  the  effect  of  kinetic  scheme  on  NO  mechanism  is 
analyzed. 


Research  methods 

Experimental  setup 

The  experiments  are  carried  out  at  a  laboratory-scale  Adelaide 
MILD  combustion  furnace  and  details  of  which  have  been 
given  in  ref.  [45,46].  The  combustion  chamber  is  well  insulated 
and  allows  only  -20%  of  the  total  heat  input  to  be  conducted 
through  the  walls.  As  a  result,  it  usually  takes  a  warm-up  time 
of  about  1.5  h  from  a  cold  state  to  steady-state  operation. 
There  are  five  openings  accommodating  interchangeable 
insulating  window  plugs  or  UV  grade  fused  silica  windows. 
The  five  openings  are  equally  spaced  vertically  down  three 
sides  of  the  furnace.  To  control  the  heat  load,  two  U-shaped 
cooling  tubes  with  variable  heat  exchange  areas  are  used.  The 
heat  exchangers  can  be  inserted  through  any  of  the  window 
openings.  For  this  investigation  their  exposed  surface  areas 
are  0.03  m2  each.  These  heat  exchangers  remove  4.46  kW  and 
4.93  kW  of  heat  on  average  for  the  firing  rates  of  10  kW  and 
15  kW,  respectively. 

The  time-averaged  furnace  temperature  is  measured  with 
a  bare,  fine-wire,  type  R  (Pt— Pt— 13%Rh)  thermocouple  of 
254  nm  diameter  wire  with  a  bead  diameter  of  1.2  mm  under 
steady-state  conditions.  Five  measurement  points  along  the 
axial  centerline  of  the  furnace  are  selected  as  temperature 
measurement  locations.  The  exhaust  temperature  is 
measured  with  a  stainless  steel  sheath  type  K  (Ni— Cr)  ther¬ 
mocouple.  Exhaust  emissions  of  NO,  N02,  C02  and  02  are 
measured  using  a  TESTO  350  XL  portable  gas  analyzer.  While 
uncertainties  of  the  mean  temperatures  measured  by  thermal 
couples  are  about  ±4  K,  the  analyzer  measurement  accuracies 
are  estimated  to  be:  [NO]  =  ±5  ppm  and  [N02]  =  ±5  ppm.  The 
analyzer  is  checked  with  a  calibration  gas  to  yield  total  NOx 
emissions  accuracies  better  than  ±5  ppm.  Total  NOx  emission 
(NO  +  N02)  concentrations  are  reported  by  volume  on  a  dry 
basis  corrected  to  3%  02  concentration.  By  correcting  to  a 
specific  02  level,  true  comparisons  of  emissions  levels  can  be 
made  because  the  effect  of  various  degrees  of  dilutions  has 
been  removed,  while  still  retaining  a  familiar  mole-fraction¬ 
like  variable. 

The  burner  consists  of  a  single  fuel  nozzle  on  the  axis  of  the 
furnace,  and  four  exhaust  and  four  air  ports  arrange  sym¬ 
metrically  in  a  ring  pattern  on  the  same  wall.  The  schematic 
figure  of  the  MILD  combustion  furnace  can  be  found  in  ref. 
[45].  The  inner  diameter  of  the  central  fuel  nozzle  and  four  air 
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nozzles  is  Df  =  7.2  and  Da  =  4.0  mm,  respectively.  In  the  pre¬ 
sent  experiments,  the  room-temperature  reactants  (  —  288  K) 
are  injected  into  the  furnace  and  the  equivalence  ratio  (<£>)  of 
0.8  is  used  and  kept  constant.  The  furnace  is  operated  with  a 
thermal  input  of  13  kW  when  only  natural  gas  (NG)  is  used  as 
fuel.  The  NG  supply  is  kept  constant.  With  the  addition  of 
hydrogen,  the  mass  fraction  of  hydrogen  in  the  fuel  jet  (NG/H2) 
is  increased  from  0%  to  14.4%  (by  mass)  and  the  thermal  input 
is  increased  from  13  kW  to  19.1  kW.  The  properties  and 
composition  of  natural  gas  used  in  the  present  experiment  is 
shown  in  Table  2. 

To  our  best  knowledge,  there  are  four  kinds  of  definitions 
of  the  MILD  combustion,  i.e.,  (1)  the  definition  of  Wiinning  and 
Wunning  [11],  (2)  the  definition  of  Kumar  et  al.  [47  ,  (3)  the 
definition  of  Cavaliere  et  al.  [8]  and  (4)  the  definition  of  the 
present  authors  [45,46,48].  Wunning  and  Wunning 
defined  the  MILD  combustion  regime  using  the  global  pa¬ 
rameters  of  gas  recirculation  rate  (Kv)  and  furnace  tempera¬ 
ture.  The  gas  recirculation  rate  (Kj  is  the  ratio  between  the 
mass  flow  rate  of  exhaust  gases  internally  recirculated  into 
the  air  and  fuel  streams  before  reaction  and  the  total  mass 
flow  rate  injected  into  the  furnace.  The  MILD  combustion 
regime  occurs  when  Kv  >  3  and  furnace  temperature  higher 
than  1100  K  [11].  However,  in  general,  it  is  extremely  difficult 
to  measure  Kv  by  experiments,  and  thus  the  judgment  of  the 
MILD  combustion  is  not  straight-forward.  Kumar  et  al.  [47] 
defined  the  MILD  combustion  based  on  the  mean  square 
temperature  fluctuation  over  the  combustor  (T'2).  T'2  is  calcu¬ 
lated  as: 


and  the  MILD  combustion  is  defined  as  (T,2/Tmean)  <  15%. 
Nevertheless,  it  is  also  not  easy  to  obtain  T'2  and  Tmea n  through 
experiments,  because  the  estimation  of  Tmean  needs  the 
temperature  data  everywhere  inside  the  entire  furnace.  The 
MILD  combustion  defined  by  Cavaliere  et  al.  [8]  is  “a  com¬ 
bustion  process  when  the  inlet  temperature  of  the  reactant 
mixture  is  higher  than  mixture  self-ignition  temperature 
whereas  the  maximum  allowable  temperature  increase  with 
respect  to  inlet  temperature  during  combustion  is  lower  than 
mixture  self-ignition  temperature.”  This  definition  is  based  on 
the  well-stirred  reactor.  The  definitions  of  the  MILD  combus¬ 
tion  of  burning  gaseous,  liquid  and  solid  fuels  are  discussed  in 
our  previous  study  48].  For  the  present  study  of  firing  gaseous 
fuels,  the  MILD  combustion  in  experiments  is  defined  as  the 
flameless  combustion  where  no  visible  flame  front  occurs  at 
all.  Although  when  the  MILD  combustion  occurs,  the  charac¬ 
teristics  of  the  thermal  and  gas  species  fields  are  different 
from  that  of  the  conventional  combustion.  It  is  better  to  use  a 


straight-forward  method  to  judge  the  occurrence  of  MILD 
combustion  in  experiments.  In  the  present  experiments, 
when  there  is  no  visible  flame,  the  measured  temperature 
fields  are  semi-uniform,  and  the  NOx  and  CO  emissions  are 
lower  than  40  ppm  and  10  ppm,  respectively.  Therefore,  for 
gaseous  fuels,  the  MILD  combustion  can  be  indeed  defined 
based  on  the  appearance  as  the  flameless  combustion  where 
no  visible  flame  front  occurs  at  all.  The  occurrence  of  the 
flameless  combustion  can  ensure  the  low  NOx  emission 
results. 


CFD  calculations 


Computational  fluid  dynamic  (CFD)  simulations  are  per¬ 
formed  to  investigate  the  NO  mechanism  of  MILD  combustion 
in  Cases  1—4.  The  modeling  detail  and  validation  can  be  found 
in  Refs.  [46,49—52]  and  only  a  brief  description  is  provided 
below. 

The  commercial  software  FLUENT  code  [53]  is  used  to 
model  these  furnace  utilizing  steady  state,  implicit,  finite 
volume  based  models.  The  pressure  velocity  coupling  is  the 
SIMPLE  algorithm.  Full  hexahedral  grid  (about  2000,000  hex- 
ahedral  cells)  is  adopted  to  minimize  the  grid  size  and 
appropriate  refinement  of  grid  is  performed  in  the  regions 
with  higher  gradients.  The  adequacy  of  this  mesh  (i.e.,  the  grid 
independence)  is  verified  by  comparing  the  results  with  those 
obtained  using  a  finer  grid  with  4,000,000  cells.  The  compari¬ 
son  between  these  two  cases  shows  a  high  consistency  in  the 
results.  The  realizable  k— e  model  with  the  standard  wall 
function  is  taken  for  modeling  the  turbulent  flow. 

The  eddy  dissipation  concept  (EDC)  [54]  is  used  with  the 
detailed  reaction  mechanism  of  GRI-Mech  2.11  [55]  to  simulate 
the  oxidation  of  the  fuel.  The  EDC  model  is  an  extension  of  the 
eddy  dissipation  model  to  consider  detailed  chemical  mech¬ 
anisms  in  turbulent  flows  and  can  capture  finite  rate  chem¬ 
istry  effects.  The  EDC  model  assumes  that  reactions  occur  at 
fine  scales.  The  reaction  at  fine  scales  is  assumed  to  occur  as  a 
constant  pressure  reactor.  The  characteristic  length  fraction 
of  fine  scales  (f)  and  the  chemical  residence  time  scale  (r)  of 
fluid  in  the  fine  structures  is  expressed  by. 


£ 


/Ue\ a/4 

\w)  ’ 


U\  1/2 
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where  Q  is  the  volume  fraction  constant  (=2.1377)  and  CT  is 
the  time  scale  constant  (=0.4083).  The  evolutions  of  species 
concentrations  are  then  computed  by  integrating  the  chem¬ 
istry  within  those  fine  scales.  In  EDC  model,  the  species  con¬ 
servation  equation  takes  the  following  general  form: 


-v;  i  +  Ri 


Table  2  -  Properties  and  composition  of  natural  gas  used  in  the  present  experiment. 

Fuel  LHVb  CH4  C2Hs 

C02  n2 

c3h8 

C4H10 

C5Hi2 

C6Hi4 

Natural  gasa  51.154  91.36  4.364 

2.084  1.278 

0.62 

0.20 

0.055 

0.04 

a  Gas  analysis  provided  by  Origin  Energy  Australia. 
b  Lower  heating  value  (MJ/kg). 
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where  Yit  J  {  and  Rj  are  respectively  the  local  mass  fraction  of 
each  species  i,  the  diffusion  flux  and  the  net  rate  of  production 
by  chemical  reaction.  The  chemistry  source  term  (R*)  is 
computed  by 

Ri  =  J(T^)(Y'-Yi)  (4) 

JfC 

with  Yj  being  the  fine-scale  species  mass  fraction  after  react¬ 
ing  over  the  time  (r).  The  evolution  of  Y*  depends  also  on  the 
chemical  kinetic  mechanism.  Note  that  the  expression  (4)  is 
only  valid  when  the  turbulence  Reynolds  number  Ret  >  65  [56]. 
If  Ret  <  65,  the  early  ignition  problem  may  occur.  This  problem 
can  be  avoided  by  modifying  the  EDC  model  constants  [56,57], 

Although  discrepancies  between  the  EDC  simulations  and 
experiments  are  notable  at  low  oxygen  levels  [42,58],  the 
similar  may  occur  for  extremely  low  oxygen  cases  if  the  time 
scale  constant  of  the  EDC  model  is  set  as  Cr  =  0.4083  [56,57] . 
Since  the  gradients  of  temperature  and  species  are  signifi¬ 
cantly  lower  in  MILD  than  conventional  combustion,  the 
former  should  have  a  greater  residence  time  in  the  fine 
structures  of  the  EDC  model  than  does  the  latter  [56,57], 
Moreover,  the  reaction  rate  of  MILD  combustion  is  far  lower 
than  that  of  conventional  combustion.  Therefore,  CT  should  be 
taken  greater  than  0.4083  so  that  the  simulated  reaction  rate  is 
reduced.  It  follows  that  the  characteristics  of  the  MILD  reac¬ 
tion  operating  at  extremely  low  oxygen  levels  can  be  captured. 
Accordingly,  in  some  of  our  modeling  studies  of  the  MILD 
combustion,  especially  under  extremely  low  oxygen  levels, 
the  time  scale  constant  is  increased  from  CT  =  0.4083  to 
Cr  =  1.5  or  3.0  [56,57]. 

For  the  present  MILD  combustion  occurring  in  a  furnace, 
since  Ret  >  65  (see  Ref.  [46]  for  detail),  the  default  values  of  EDC 
model  constants  of  Q  =  2.1377  and  Cr  =  0.4083  are  adopted.  As 
shown  in  Section  3.2,  the  predicting  temperature  distributions 
agree  well  with  experiments. 

Differential  diffusion  is  considered  by  representing  mo¬ 
lecular  diffusion  coefficients  for  each  species  as  a  fourth-order 
polynomial  function  of  temperature.  The  discrete  ordinate 
radiation  model  is  adopted  with  a  weight  sum  of  gray  gas 
model  to  model  the  radiation.  Once  the  flow  and  thermal 
fields  have  been  obtained  from  the  model  and  validated  with 
the  experiment,  the  NO  route  paths  are  calculated  as  a  post¬ 
processing  operation  in  FLUENT  [53].  The  NO  production 
typically  appears  in  low  concentrations  and  the  NO  chemistry 
has  a  negligible  influence  on  the  predicted  flow  field,  tem¬ 
peratures  and  major  combustion  product  concentrations  and 
thus  the  post-processing  operation  of  the  NO  calculation  is 
reasonable.  The  thermal,  prompt,  N20-intermediate,  NNH  and 
NO-reburning  mechanisms  are  considered  in  the  present 
simulation.  The  O  radical  concentration  in  thermal  NO  and 
N20-intermediate  routes,  as  well  as  the  OH  radical  concen¬ 
tration  in  thermal  NO,  are  predicted  using  the  EDC  model  with 
GRI-Mech  2.11  mechanism  (this  method  is  termed  as  the 
instantaneous  approach  in  FLUENT  53]).  The  prompt  NO 
formation  is  modeled  following  De  Soete  [59].  The  N20-inter- 
mediate  mechanism  is  assumed  at  the  quasi-steady-state  53]. 
The  thermal,  prompt  and  N20-intermediate  routes  are  based 
on  kinetic  mechanisms  with  Arrhenius  equations  integrated 
over  a  probability  density  function  (PDF)  for  temperature  in 


order  to  take  into  account  the  effect  of  turbulent  fluctuations 
on  the  mean  reaction  rates.  For  the  NNH  route,  because  there 
is  no  commercial  code  available,  we  use  the  EDC  model  with 
the  NNH  mechanism  [60]  to  simulate  its  production,  see  ref. 
[45,46]  for  details. 

Boundary  conditions  are  velocity  inlet,  pressure  outlet  and 
temperature  wall  boundary.  The  temperature  of  wall  bound¬ 
ary  is  set  by  using  the  user  defined  function  of  FLUENT  53]  to 
be  consistent  with  the  experimental  measurement.  The 
second-order  scheme  is  employed  for  pressure  and  the 
second-order  upwind  scheme  for  momentum,  turbulent  ki¬ 
netic  energy,  turbulent  dissipation,  species  transport  and  en¬ 
ergy  equations.  Convergence  is  obtained  when  residuals  are 
below  1CT6  for  the  energy  and  1CT5  for  all  other  variables.  The 
outlet  temperature  and  velocity  are  monitored  and  their  var¬ 
iations  are  allowed  to  be  within  1  K  and  0.1  m/s,  respectively, 
to  achieve  the  convergence  of  their  solution. 

Chemical  kinetics  calculations 

MILD  combustion  is  characterized  by  intense  internal  recir¬ 
culation  and  low  temperature  gradient  inside  the  furnace. 
Theoretically,  the  entire  in-furnace  region  of  the  MILD  com¬ 
bustion  furnace  can  be  simplified  as  a  large  well-stirred 
reactor  (WSR),  because  chemical  timescales  are  much  larger 
than  turbulent  timescales  due  to  the  relatively  strong  internal 
recirculation  and  dilution  of  the  reactant  [24,61].  De  Joannon 
et  al.  [61]  reported  the  similarity  between  concepts  of  MILD 
condition  and  WSR.  They  used  the  WSR  to  investigate  the 
characteristics  of  MILD  combustion  regime.  Mardani  et  al.  [24] 
also  adopted  the  WSR  to  explore  the  CO  and  C02  formation  in 
CHVH2  blended  combustion  under  MILD  condition.  In  addi¬ 
tion,  the  EDC  model  assumes  that  reaction  occurs  in  small 
turbulent  structures  (called  the  fine  scales)  and  the  fine  scale 
are  assumed  as  isobaric,  adiabatic  and  perfectly  stirred  re¬ 
actors  [54],  Recently,  Minamoto  and  Swaminathan 
discovered  through  direct  numerical  simulation  that  the 
WSR-based  model  is  appropriate  for  the  modeling  of  MILD 
condition.  Therefore,  there  are  some  qualitative  similarities 
between  characteristics  of  the  MILD  regime  and  the  WSR.  For 
the  present  study,  the  chemical  kinetics  calculation  is  carried 
out  with  the  WSR  model  to  investigate  the  NO  mechanism. 
The  WSR  analysis  is  performed  by  solving  the  species  con¬ 
servations  equations  under  constant  reactor  temperature 
over  a  residence  time  using  CHEMKIN  4.1  [63]. 

The  residence  time  (r)  is  set  as  1.00  s  in  the  WSR  model.  In 
our  preliminary  calculation,  the  influence  of  r  on  the  reaction 
pathways  of  NO  mechanism  is  examined.  It  is  noticed  that  r 
has  significant  influence  on  the  elementary  reaction  path¬ 
ways,  when  r  is  varied  from  0.01  s  to  1.00  s.  The  reason  is  that 
the  NO  formation  reaction  is  slow.  By  setting  r  =  1.00  s  it  is 
possible  to  ensure  that  the  NO  reaction  path  is  near  the 
chemical  equilibrium  state.  It  is  worth  noting  that  other  re¬ 
searchers  [64,65]  also  adopted  the  residence  time  in  the  range 
from  0.30  s  to  1.00  s  to  investigate  the  NO  formation  and 
reburning  chemistry. 

It  is  vital  to  choose  an  appropriate  chemical  mechanism  for 
the  WSR  calculation.  GRI-Mech  2.11  is  suitable  for  modeling 
the  condition  of  roughly  1000—2500  K,  10  torr  to  10  atm,  and 
equivalence  ratio  from  0.1  to  5  [55].  Previous  studies  [66-69] 
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showed  that  the  full  chemical  mechanism  of  GRI-Mech  2.11  is 
better  than  GRI-Mech  3.0  in  the  modeling  of  the  NO  mecha¬ 
nism.  Therefore,  GRI-Mech  2.11  is  adopted  in  the  present 
modeling. 

The  sensitivity  analysis  is  performed  to  clarify  the  main 
reaction  pathways  for  NO  formation  and  reduction.  The 
sensitivity  analysis  provides  the  values  of  the  first-order 
sensitivity  coefficients  (A/Yj)  (<5Y;-/<5Ai).  Here,  A*  is  the  pre¬ 
exponential  factor  for  reaction  i,  and  Yj  is  the  mass  fraction 
for  the  jth  species.  Thus,  the  sensitivity  coefficients  can  be 
interpreted  as  the  relative  change  in  the  predicted  concen¬ 
tration  for  species;  caused  by  increasing  the  rate  constant  for 
reaction  i. 

The  MILD  cases  investigated 

Table  3  contains  a  summary  of  all  MILD  condition  cases.  Cases 
1—4  are  investigated  experimentally  and  numerically.  For 
these  cases,  the  furnace  temperature,  extracted  heat  and 
exhaust  NOx  emissions  are  measured  in  the  experiments, 
while  the  in-furnace  flow  field  and  NO  emission  are  obtained 
by  CFD.  The  effects  of  hydrogen  addition  (Cases  5-7),  reactor 
temperature  (Cases  8-10),  oxygen  level  (Case  11-13)  and 
equivalence  ratio  (Cases  14—16)  on  the  NO  formation  and 
reduction  mechanisms  are  explored  by  chemical  kinetics 
calculations  with  the  well-stirred  reactor  (WSR)  model.  For 
Cases  11—13,  the  initial  oxygen  mass  fraction  is  changed  and 
the  equivalence  ratio  is  kept  constant  by  varying  the  N2 
diluting  level.  For  Case  16,  the  fuel  rich  case  is  examined 
because  there  is  fuel  rich  region  existing  inside  the  practical 
furnace. 


Results  and  discussion 

Experimental  furnace  results 

In  Fig.  1  the  effects  of  adding  hydrogen  to  the  natural  gas  fuel 
stream  in  the  MILD  combustion  furnace  are  shown.  The 


hydrogen  addition  varied  from  0  to  14.4%  as  follows;  (a) 
Yh2  =  0%,  (b)  y{j2  =  5.7%  and  (c)  y£2  =  14.4%.  The  thermal  field 
is  obtained  by  measuring  the  mean  temperature  using  a 
thermocouple,  while  the  simplified  flow  field  is  deduced  from 
the  CFD  results  and  represented  by  the  dashed  lines  on  the 
left. 

For  the  three  cases  of  Fig.  1,  the  in-furnace  temperature  is 
semi-uniform  and  relatively  low  (<1300  K).  The  temperature 
difference  across  the  furnace  is  small,  mainly  less  than  50  K  in 
each  case.  Moreover,  as  the  hydrogen  is  added  to  natural  gas, 
although  the  furnace  temperature  is  increased,  the  location  of 
the  high  temperature  zone  does  not  change.  From  Fig.  1,  it  is 
clear  that  the  high  temperature  region  is  located  at  the  junc¬ 
tion  of  the  upward  fresh  fuel  jet  and  downward  recirculated 
air  jet.  The  flow  rate  of  the  air  jet  is  ten  times  greater  than  the 
fuel  jet  and  thus  the  jet  momentum  of  the  former  is  signifi¬ 
cantly  higher  than  the  latter.  Therefore,  the  strong  air  jet 
dominates  the  flow  and  entrains  the  weak  fuel  jet  at  the  lower 
part  of  the  furnace.  There  is  a  V-shaped  junction  region  of  the 
fuel  and  air  jets.  The  MILD  combustion  mainly  occurs  in  this 
V-shaped  region  and  thus  the  temperature  is  relatively  high  in 
this  zone  (see  Fig.  1).  The  temperature  in  the  main  reaction 
region  increases  with  the  addition  of  hydrogen.  For  instance, 
at  y{j2  =  0%,  the  temperature  at  the  main  reaction  zone  is 
approximately  1200  K;  at  y£  =  5.7%,  the  temperature  at  the  V- 
shaped  zone  is  1250  K;  at  Y^  =  14.4%,  it  is  1300  K.  Note  that, 
although  ig.  1  only  shows  the  two-dimensional  cross  section 
of  the  furnace,  this  section  incorporates  the  reactants  jets  and 
so  represents  the  most  important  reaction  zone  in  the 
furnace. 

It  is  also  worth  noting  that  the  temperature  distribution  of 
the  MILD  combustion  in  the  present  experiment  is  different 
from  our  previous  study  using  the  same  furnace,  see  ref.  [45] 
for  details,  because  the  burner  configuration  is  different.  For 
our  previous  study,  the  burner  consists  of  a  central  air  jet  and 
four  side  fuel  jets.  Therefore,  the  junction  of  the  previous  fuel 
and  air  jets  was  located  around  the  furnace  central  region  and 
thus  high  temperature  is  also  over  there.  For  the  present 
study,  a  V-shaped  junction  region  is  formed  by  the  side  air  jets 


Table  3  -  Summary  of  the  MILD  cases  investigated. 

Cases 

Fuel 

YfH2  (%) 

Yg2  (%) 

0 

EXP  &  CFD 

WSRa 

f  WSR 

(Yh2,  YCh4,  Yq2,  Yn2)  in  WSR 

l 

NG 

0.0 

23.0 

0.8 

V 

— 

— 

— 

2 

NG  +  H2 

5.7 

23.0 

0.8 

V 

— 

— 

— 

3 

NG  +  H2 

10.8 

23.0 

0.8 

V 

— 

4 

NG  +  H2 

14.4 

23.0 

0.8 

V 

— 

5 

ch4 

0.0 

3.0 

0.8 

— 

V 

1345 

(0.000%,  0.597%,  2.982%,  96.421%) 

6 

ch4  +  h2 

10.0 

3.0 

0.8 

— 

V 

1345 

(0.054%,  0.488%,  2.984%,  96.474%) 

7 

ch4  +  h2 

20.0 

3.0 

0.8 

— 

V 

1345 

(0.100%,  0.398%,  2.985%,  96.517%) 

8 

ch4  +  h2 

20.0 

3.0 

0.8 

— 

V 

1100 

(0.054%,  0.488%,  2.984%,  96.474%) 

9 

ch4  +  h2 

20.0 

3.0 

0.8 

— 

V 

1345 

(0.054%,  0.488%,  2.984%,  96.474%) 

10 

ch4  +  h2 

20.0 

3.0 

0.8 

— 

V 

1600 

(0.054%,  0.488%,  2.984%,  96.474%) 

11 

ch4  +  h2 

20.0 

3.0 

0.8 

— 

V 

1345 

(0.054%,  0.488%,  2.984%,  96.474%) 

12 

ch4  +  h2 

20.0 

6.0 

0.8 

— 

V 

1345 

(0.108%,  0.971%,  5.935%,  92.986%) 

13 

ch4  +  h2 

20.0 

9.0 

0.8 

— 

V 

1345 

(0.161%,  1.229%,  8.855%,  89.535%) 

14 

ch4  +  h2 

20.0 

3.0 

0.5 

— 

V 

1345 

(0.034%,  0.306%,  2.990%,  96.670%) 

15 

ch4  +  h2 

20.0 

3.0 

1.0 

— 

V 

1345 

(0.068%,  0.609%,  2.980%,  96.343%) 

16 

ch4  +  h2 

20.0 

3.0 

1.5 

— 

V 

1345 

(0.101%,  0.911%,  2.970%,  96.018%) 

a  r  and  p  for  WSR  is  1  s  and  1  atm,  respectively. 
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Fig.  1  —  Measured  furnace  temperature  distribution  under  the  MILD  combustion  at  x  =  0,  together  with  a  simplified  diagram 
of  the  flow  field  obtain  from  CFD.  (a)  y£2  =  0%,  (b)  y£2  =  5.7%  and  (c)  y{j2  =  14.4%  (Cases  1,  2,  and  4). 


and  central  fuel  jet  configuration.  Hence  the  V-shaped  high 
temperature  region  is  formed  there  as  shown  in  Fig.  1. 

Fig.  2a  shows  the  effects  of  hydrogen  addition  on  the 
exhaust  mass  fraction  of  02  (Y0z)  and  C02  (YCo2),  as  well  as  the 
NOx  emission  (XNOx,  ppm  in  volume)  while  Fig.  2b  displays  the 
furnace  average  temperature  (Taverage)>  exhaust  temperature 
(Texhaust),  thermal  input  (Qinput)  and  extracted  heat  (Qextracted) 
versus  the  hydrogen  addition.  As  the  mass  fraction  of 
hydrogen  in  the  fuel  (y{j  )  is  increased  from  0  to  14.4%,  Qinput 
increases  from  13.0  kW  to  19.1  kW,  Qextracted  increases  from 
5.3  kW  to  6.3  kW,  and  Taverage  increases  from  1164  K  to  1232  K. 


Note  that  Taverage  is  obtained  by  averaging  temperatures  along 
the  furnace  centreline.  As  y{j  increases  from  0  to  5.7%,  the 
NOx  emission  increases  from  18  ppm  to  35  ppm.  However, 
interestingly,  as  Y^  is  further  increased  from  5.7%  to  14.4%, 
the  NOx  emission  is  not  changed  and  remains  constant  at 
35  ppm.  Namely,  as  y£  is  increased  from  5.7  to  14.4%,  both  the 
Taverage  and  Qextracted  are  increased,  but  the  NOx  emission  re¬ 
mains  unchanged.  This  cannot  be  explained  by  the  conven¬ 
tional  thermal  NO  mechanism  because  the  previous 
knowledge  tells  that  the  NOx  emission  is  generally  increased 
as  the  furnace  temperature  increases  [10,70].  Therefore,  it  is 


Fig.  2  -  Effects  of  hydrogen  addition  on  (a)  the  exhaust  emissions  of  02,  C02,  and  NOx  and  (b)  the  furnace  average 
temperature,  exhaust  temperature,  thermal  input,  and  extracted  heat  (Cases  1-4). 
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necessary  to  investigate  the  NO  formation  mechanism  under 
the  MILD  combustion  of  hydrogen-containing  hybrid  fuels  in 
detail,  which  is  presented  in  next  section.  It  is  worth 
mentioning  that  the  CO  emissions  for  all  experimental  cases 
are  too  small,  below  the  threshold  of  the  gas  analyzer,  and 
thus  not  recorded.  The  extremely  low  emissions  of  CO  for  all 
experimental  cases  indicate  that  the  fuel  is  completely 
oxidized  in  the  experimental  combustion,  regardless  of  the 
hydrogen  addition. 


r^H2  0 

/m  =  5.7% 

/h2=  10.8% 

=  14.4% 


CFD  results 

In  this  Section,  CFD  results  are  presented  in  order  to  inves¬ 
tigate  the  mechanisms  of  NO  formation  and  destruction  due 
to  the  addition  of  hydrogen.  It  is  well  established  that  the 
temperature  is  vital  for  the  evolution  of  NO,  especially  that 
some  NO  mechanisms  are  strongly  dependent  on  tempera¬ 
ture.  For  the  present  study,  the  in-furnace  and  wall  tem¬ 
peratures  of  the  CFD  and  experiment  are  compared  firstly, 
and  then  the  NO  mechanism  calculation  is  carried  out.  In 
Fig.  3a— c  the  comparison  of  the  experimental  and  numerical 
in-furnace  temperature,  obtained  along  the  z  axis  at 
x  =  0  mm  and  y  =  0  mm,  50  mm  and  100  mm  of  Case  2  is 
shown.  It  is  clear  that  the  CFD  results  agree  well  with  the 
measured  average  temperatures.  Fig.  3d  shows  a  comparison 
of  the  experimental  and  numerical  furnace  wall  tempera¬ 
ture,  obtained  along  the  z  axis  at  x  =  0  mm  and  y  =  140  mm. 
The  furnace  wall  temperature  in  the  CFD  models  is  set  by 
using  the  user-defined-function  of  FLUENT  [53]  and  it  is 
consistent  with  those  measured  experimentally.  Therefore, 
the  thermal  field  of  the  CFD  result  is  validated  with  the  ex¬ 
periments  and  provides  confidence  in  performing  analyses  of 
the  NO  field. 


Fig.  4  —  Calculated  relative  importance  of  NO  emissions  via 
thermal,  prompt,  N20 -intermediate,  and  reburning  routes 
at  the  furnace  exit  (Cases  1-4). 


Fig.  4  depicts  calculated  NO  emissions  via  thermal,  prompt, 
N20-intermediate,  and  reburning  routes  at  the  furnace  exit. 
The  calculation  method  is  introduced  in  sub-section  CFD 
calculations  .  It  is  evident  that  the  N20-intermediate  route 
dominates  the  NO  formation  in  all  the  cases,  with  its  contri¬ 
bution  growing  slightly  as  y{j  is  increased.  When  y{j  =  0%  or 
burning  natural  gas,  the  N20-intermediate  mechanism  ac¬ 
counts  for  71%  of  the  total  NO  formation.  The  rest  of  NO  is 
formed  from  the  prompt  and  NNH  routes.  The  NO-reburning 
reaction  is  strong  and  reduces  the  NO  formation  by  about 
30%.  This  is  consistent  with  those  of  Li  et  al.  [45,46]  who  used 
different  burner  configurations  and  fuel  mixture  from  the 
present.  Those  authors  found  that  the  NO  mainly  comes  from 
the  N20-intermediate  route,  regardless  of  the  hydrocarbon 
fuel  type  and  oxidant  type  [45]  and  premixing  pattern  [46]. 
Therefore,  it  is  concluded  that,  when  0  <  0.8,  the  N20-inter- 
mediate  route  at  these  temperatures  is  always  the  dominant 
NO  formation  mechanism  under  MILD  combustion,  regardless 
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• 

Experiments 

CFD 

(b)  x  =  0  mm,  y  =  50  mm 

J _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ u 


100  200  300  400  500  600  100  200  300  400  500  600 


z  [mm]  z  [mm] 

Fig.  3  —  Comparison  of  the  experimental  and  numerical  furnace  temperature  profiles  along  the  z  axis  at  x  =  0  mm  and 
y  =  0  mm,  50  mm,  100  mm  and  140  mm  (Case  2). 
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of  the  hydrocarbon  fuel  type,  the  hydrogen  addition,  the 
oxidant  type,  and  burner  configuration. 

As  y{j2  is  increased  from  0  to  5.7%,  the  N20-intermediate 
yields  more  NO  formation  and  the  NO  reduction  from  NO- 
reburning  is  reduced.  So,  the  measured  NO  emission  is 
increased  (see  rig.  2).  However,  as  y{j  is  further  increased 
from  5.7%  to  14.4%,  NO  formation  from  the  prompt  route  is 
decreased  while  that  from  the  NNH  route  is  increased.  As  a 
result,  the  total  NO  emission  changes  insignificantly  as 
measured  in  experiments  (  lg.  2).  Fig.  4  also  reveals  that  un¬ 
surprisingly  at  these  temperatures,  the  thermal  NO  mecha¬ 
nism  can  be  ignored  for  all  the  MILD  cases,  which  differs 
significantly  from  the  conventional  combustion.  To  more 
deeply  understand  the  NO  mechanism  under  the  MILD  com¬ 
bustion  of  CH4/H2,  detailed  chemical  kinetic  calculations 
using  WSR  model  are  carried  out  to  examine  the  effects  of 
hydrogen  addition,  reactor  temperature,  mass  fraction  of  ox¬ 
ygen  and  equivalence  ratio,  and  the  results  are  presented 
below. 

Chemical  kinetics  calculation  results 

Effects  ofH2  addition 

Fig.  5  displays  the  reaction  path  diagram  for  reactive  nitrogen 
conversion  in  the  MILD  combustion  of  CH4  and  CH4/H2  mix¬ 
tures.  N2  convert  to  NO  mainly  through  the  N2  — ►  N20  — ►  NO, 
N2  NNH  ->  NO,  N2  HCN  ->  NCO  NO  and  N2  -►  N  -►  NO 
routes.  The  N2  — ►  N20  — ►  NO  pathway  corresponds  to  the  N20- 
intermediate  mechanism.  The  N2  — ►  NNH  — ►  NO  trail  is  the 
NNH  mechanism  while  the  prompt  NO  route  is 
N2  — ►  HCN  — ►  NCO  —>  NO.  Therefore,  it  is  consistent  with  the 
CFD  results  of  the  above  sub-section,  the  reaction  path  dia¬ 
gram  analysis  also  shows  that  N20-intermediate,  NNH  and 


prompt  routes  are  the  vital  reactions  in  producing  NO  under 
the  MILD  condition.  Moreover,  the  NO  is  reduced  through  the 
NO  ->  HNO,  NO  HCNO,  NO  HNCO  and  NO  HCN 
pathways.  This  indicates  that  the  NO-reburning  is  important 
under  the  MILD  condition,  reducing  the  NO  formation  mainly 
through  reactions  CH2  +  NO  H  +  HNCO, 

H  +  NO  +  M  HNO  +  M,  CH2  +  NO  ^  H  +  HCNO,  and 
CH2  +  NO  OH  +  HCN.  Furthermore,  the  mutual  conversion 
between  NO  and  N02  is  significant.  The  N0/N02  conversion 
reactions  become  strong  in  low  temperatures  71].  However, 
for  conventional  combustion,  no  N02  is  formed  in  high- 
temperature  regions  because  N02  destruction  reactions  are 
extremely  strong  at  high  temperatures  71]. 

The  results  for  Cases  5—7  of  the  first-order  sensitivity 
analysis  for  NO  are  presented  in  Fig.  6.  The  first  two  sensitivity 
reactions  for  the  NO  reduction  are  N20  +  H  ^  N2  +  OH  and 
H  +  02  +  N2  H02  +  N2.  The  sensitivity  coefficients  for  these 
reactions  are  negative,  indicating  that  increasing  the  reaction 
rate  will  lead  to  more  NO  reduction.  H02  is  largely  formed 
under  the  low  temperature  regime  of  MILD  condition  through 
this  reaction  and  thus  significant  NO  is  converted  to  N02  via 
reaction  NO  +  H02  N02  +  OH.  The  first  three  maximum 

sensitivity  reactions  for  the  NO  formation  are 
N20  +  H  NH  +  NO,  N20  +  O  2NO  and  N20  (+M)  N2  +  O 
(+M).  The  sensitivity  coefficients  of  these  reactions  are  posi¬ 
tive  which  indicates  that  increasing  the  rate  of  these  reactions 
will  lead  to  more  NO  production.  Reactions 
N20  +  H  NH  +  NO  and  N20  +  O  2NO  correspond  to  the 
N20-intermediate  route.  Therefore,  the  N20-intermediate  re¬ 
action  is  the  most  important  route  for  the  NO  production 
under  these  MILD  condition.  The  NNH  route 
(NNH  +  O  ^  NH  +  NO)  is  significant  as  well  but  it  is  less 
important  than  the  N20-intermediate  route.  Fig.  6  further 


Fig.  5  -  Reaction  path  diagram  for  reactive  nitrogen  conversion  in  the  MILD  combustion  of  CH4  and  CH4/H2  mixtures.  The 
dashed  lines  denote  pathways  only  important  in  the  MILD  combustion  of  CH4  (Case  5).  The  red  lines  denote  those  pathways 
becoming  more  important  in  the  MILD  combustion  of  CH4/H2  mixtures  (Cases  6-7).  (For  interpretation  of  the  references  to 
color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  6  -  Results  of  the  first-order  NO  sensitivity  analysis  for  Cases  5-7. 


suggests  that  the  prompt  route  produces  some  NO  via  its  key 
step  CH  +  N2  HCN  +  N.  The  prompt  route  is  weaker  than  the 
N20  and  NNH  routes. 

Moreover,  as  more  hydrogen  is  added,  the  sensitivity  co¬ 
efficient  of  reaction  NNH  +  0^  NH  +  NO  is  increased  because 
more  NNH  radical  is  formed  via  reaction  N2  +  H  (+M)  NNH 
(+M).  Consequently,  thus  more  NO  is  produced  from  the  NNH 
route  with  hydrogen  addition.  However,  as  more  hydrogen 
is  added,  the  sensitivity  coefficients  of  reaction 


CH  +  N2  HCN  +  N  is  decreased  due  to  the  consumption  of 
CH  through  reaction  CH  +  H2  ++  H  +  CH2  Therefore,  NO 
formed  from  the  prompt  route  is  decreased  with  the  hydrogen 
addition.  This  WSR  result  is  consistent  with  the  CFD  result  of 

Fig.  4. 

Effects  of  reactor  temperature 

Fig.  displays  the  first-order  sensitivity  analysis  for  NO  for 
Cases  8-10.  Clearly,  the  N20  route  (N20  +  H  ^  NH  +  NO  and 
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Fig.  7  -  First-order  sensitivity  analysis  for  NO  for  Cases  8-10. 
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N20  +  O  2NO)  is  dominant,  regardless  of  the  reactor  tem¬ 
perature.  The  NNH  (NNH  +  O  NH  +  NO)  and  prompt 
(CH  +  N2  ^  HCN  +  N)  routes  are  less  important  but  also  play 
their  role  in  production  NO.  The  NO-reburning  reaction 
CH2  +  NO  H  +  HNCO  reduces  NO  production.  As  the  reactor 
temperature  is  increased  from  1100  K  to  1600  K,  the  impor¬ 
tance  of  the  NNH  route  decreases,  while  those  of  prompt  and 
thermal  routes  increase.  Importantly,  as  TWSR  is  reduced  from 
1600  K  to  1100  K,  the  sensitivity  of  reactions 
H  +  02  +  N2  H02  +  N2  and  H  +  02  +  H20  H02  +  H20 
decreases.  H02  is  formed  through  these  reactions  and  thus 
more  NO  is  converted  to  N02  through  reaction 
NO  +  H02  N02  +  OH.  Consequently,  the  relative  importance 

of  N02  to  the  total  NOx  is  increased  under  low  temperature 
condition  found  in  MILD  condition. 

Effects  of  the  mass  fraction  of  oxygen 

Fig.  8  illustrates  the  first-order  sensitivity  analysis  for  NO  for 
Cases  11-13.  Clearly,  as  the  mass  fraction  of  oxygen  (Yg  )  is 
increased,  the  importance  of  the  N20  route  increases  because 
more  NO  is  produced  from  the  reaction  N20  +  O  2NO. 
However,  the  importance  of  the  NNH  route  decreases  because 
less  NNH  radicals  are  produced  with  the  oxygen  addition  via 
reactions  N2  +  H  (+M)  NNH  (+M).  Moreover,  the  sensitivity 
of  the  prompt  mechanism  also  reduces  because  less  CH  rad¬ 
icals  are  formed  as  the  reaction  condition  becomes  fuel  lean. 

Effects  of  the  equivalence  ratio 

The  first-order  sensitivity  analysis  for  NO  for  Cases  14—16  is 
shown  in  Fig.  9.  For  0  =  0.5,  the  NO  is  mainly  produced  from 
the  N20-intermediate  route  through  reactions 
N20  +  O  2NO.  Reactions  H  +  02  +  N2  H02  +  N2  and 
H  +  02  +  H20  H02  +  H20  have  highly  sensitivity  for 


reducing  NO  via  the  reaction  NO  +  H02  N02  +  OH.  For  0  =  1, 
the  NO  is  mostly  formed  from  the  prompt  mechanism  (via  the 
reaction  CH  +  N2  HCN  +  N)  while  the  NNH  route 
(NNH  +  O  NH  +  NO)  is  also  notable.  The  NO-reburning  re¬ 
action  (CH2  +  NO  H  +  HNCO)  becomes  important  at  0  =  1. 
For  0  =  1.5,  both  the  prompt  and  NO-reburning  routes 
(CH2  +  NO  H  +  HNCO  and  HCCO  +  NO  HCNO  +  CO) 
become  extremely  stronger. 


Quantitative  analysis  of  the  NO  mechanism 

Relative  importance  of  the  NO  mechanism  through  W SR 
modeling 

The  above  analysis  on  the  NO  formation  and  destruction 
mechanisms  is  conducted  using  the  WSR  model  with  the  GRI- 
Mech  2.11  mechanism.  The  sensitivity  of  the  different  re¬ 
actions  that  contribute  to  the  NO  production  and  destruction 
is  investigated,  while  the  contribution  of  each  NO  mechanism 
to  total  NO  formation  is  not  obtained.  This  is  quantitatively 
investigated  in  this  section. 

It  is  worth  mentioning  that  reactions  of  each  NO  mecha¬ 
nism  are  complicated  and  coupled  with  each  other.  Results 
may  not  be  true  if  different  NO  pathway  are  decouple  72].  For 
the  present  investigation,  the  relative  importance  of  the 
thermal,  prompt,  N20-intermediate,  NNH  and  NO-reburning 
routes  is  obtained  by  performing  calculations  with  each 
mechanism  at  a  time.  For  each  calculation  of  a  particular 
mechanism,  relevant  reactions  are  included  to  avoid  omitting 
some  coupled  reactions.  Table  4  shows  main  reactions  of 
different  NO  routes  for  the  present  calculation. 
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Fig.  8  -  First-order  sensitivity  analysis  for  NO  for  Cases  11-13. 
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Fig.  9  —  First-order  sensitivity  analysis  for  NO  for  Cases  14-16. 


For  the  thermal  route,  the  extended  Zeldovich  mechanism 
is  considered.  For  the  prompt  route,  CH  and  CH2  radicals  react 
with  molecular  nitrogen  to  from  HCN  through  reactions 
CH  +  N2  HCN  +  N  and  CH2  +  N2  HCN  +  NH.  HCN  and  N 


convert  to  NO  via  series  of  intermediate  reactions  of  Table  4 
[73  .  For  the  N20-intermediate  route,  N2  converted  to  N20 
and  then  NO  via  reactions  N2  +  O  +  (+M)  N20  (+M), 

H  +  N20  NO  +  NH,  O  +  N20  2NO  and  other  series  of 


Table  4  -  Main  reactions  of  different  NO  routes. 

NO  formation 

Main  reactions 

Thermal 

N  +  NO  N2  +  O 

N  +  02  -  NO  +  O 

N  +  OH  ^  NO  +  H 

Prompt 

N  +  OH  ^  NO  +  H 

NH  +  H  —  N  +  H2 

CN  +  02  NCO  +  O 

NCO  +  O  ^  NO  +  CO 

NCO  +  H  ^  NH  +  CO 

HCN  +  M  —  CN  +  H  +  M 

HCN  +  O  ^  NCO  +  H 

HCN  ^  CN  +  H 

C  +  N2  ^  CN  +  N 

CH  +  N2  HCN  +  N 

CH2  +  N2  HCN  +  NH 

CH2(S)  +  N2  NH  +  HCN 

N20-intermediate 

N20  +  O  ^  N2  +  02 

N20  +  O  2NO 

N20  +  H  N2  +  OH 

N20  +  OH  N2  +  H02 

N20  (+M)  N2  +  O  (+M) 

NH  +  O  ^  NO  +  H 

NH  +  NO  N20  +  H 

NCO  +  02  ^  NO  +  C02 

NCO  +  O  NO  +  CO 

NCO  +  OH  ^  NO  +  H  +  CO 

NNH  Route 

NH  +  O  ^  NO  +  H 

NH  +  H  N  +  H2 

NH  +  OH  ^  HNO  +  H 

NH  +  OH  N  +  H20 

NH  +  02  HNO  +  O 

NH  +  02  ^  NO  +  OH 

NH  +  N  ^  N2  +  H 

NH  +  H20  ^  HNO  +  H2 

NH  +  NO  ^  N2  +  OH 
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intermediate  reactions  [65].  For  the  NNH  route,  NO  is  pro¬ 
duced  through  reactions  N2  +  H  (+M)  NNH  (+M), 
NNH  +  O  NO  +  NH  and  other  series  of  intermediate  re¬ 
actions  73].  For  the  NO-reburning  mechanism,  Dagaut  et  al. 
[74]  found  that  NO-reduction  by  the  CH4/H2  fuel  occurs  via 
HCCO  +  NO  HCN  +  C02,  HCCO  +  NO  HCNO  +  CO, 
HCNO  +  H  HCN  +  OH,  HCN  +  O  NH  +  CO, 

HCN  +  O  ^  CN  +  OH,  HCN  +  O  ^  NCO  +  H, 

HCN  +  OH  HOCN  +  H,  HCN  +  OH  ^  HNCO  +  H, 
HCN  +  OH  ^  NH2  +  CO,  CN  +  02  NCO  +  O, 

CN  +  OH  ^  NCO  +  H,  NCO  +  H  ^  NH  +  CO, 

NH  +  NO  N20  +  H  and  N20  +  H  N2  +  OH.  Other  inter¬ 
mediate  reactions  of  NO-reburning  mechanism  are  also 
considered  in  "able  4.  The  calculated  relative  importance  of 
NO  formations  and  reductions  via  thermal,  prompt,  N20-in- 
termediate,  NNH,  and  reburning  routes  of  the  WSR  model  of 
Cases  5-16  is  displayed  in  Fig.  10. 

Clearly  from  rig.  10,  for  0  <0.8  (Cases  5—14),  the  N20-in- 
termediate  route  plays  the  most  significant  role  in  the  total  NO 
formation,  regardless  of  the  hydrogen  addition  (y{j  =  0%— 
20%),  initial  mass  fraction  of  oxygen  in  oxidant  (Yg2  =  3%— 9%) 
and  reactor  temperature  (TWSR  =  1100  K-1600  K).  The  N20- 
intermediate  contributes  60%— 90%  of  the  total  NO  formation. 
The  prompt  and  NNH  routes  are  also  important.  The  NO- 
reburning  reaction  is  strong  which  reduces  30%— 50%  of  the 
total  NO  formation. 

For  Cases  5—7,  with  the  addition  of  hydrogen,  the  impor¬ 
tance  of  the  NNH  routes  increases.  However,  the  prompt  and 
NO-reburning  reactions  significantly  become  weak  with  the 
hydrogen  addition.  This  result  is  consistent  with  CFD  results 

of  Fig.  4. 

For  Cases  8-10,  effects  of  the  reactor  temperature  (TWSR)  on 
the  NO  mechanism  is  obtained.  Clearly,  the  N20-intermediate 
route  controls  the  total  NO  formation.  However,  as  TWSR  is 
increased  from  1100  K  to  1600  K,  the  relative  importance  of  the 


N20-intermediate  route  to  the  total  NO  formation  decreases 
from  93%  to  62%,  while  that  of  the  thermal  route  increases 
from  0.03%  to  24%.  For  TWSR  =  1345  K,  the  relative  importance 
of  NNH  and  NO-reburning  routes  are  higher  than  that  of 
TWSr  =  1100  K  or  1600  K. 

For  Cases  11—13,  effects  of  initial  mass  fraction  of  oxygen 
(Yq  )  on  the  NO  mechanism  is  investigated.  Obviously,  the 
N20-intermediate  mechanism  also  controls  the  total  NO  pro¬ 
duction,  regardless  of  Yg2.  As  Yg2  is  increased  from  3%  to  9%, 
the  relative  importance  of  the  N20-intermediate  increases, 
while  that  of  prompt  and  NNH  routes  decrease. 

For  Cases  14—16,  effects  of  equivalence  ratio  (0)  on  the  NO 
mechanism  is  examined.  For  Case  14,  0  =  0.5,  the  N20-inter- 
mediate  route  is  extremely  strong  and  it  accounts  for  95%  of 
the  total  NO  formation.  The  prompt  and  NNH  routes  are  un¬ 
important  at  0  =  0.5.  However,  for  Cases  15  and  16,  0  =  1  and 
1.5,  respectively.  The  prompt  route  is  so  strong  that  it  con¬ 
tributes  more  than  55%  of  the  total  NO  formation.  Especially, 
the  prompt  route  accounts  for  90%  of  the  total  NO  formation  at 
0  =  1.5.  Moreover,  the  NO-reburning  reaction  reduces 
approximately  20%  of  the  total  NO  formation  at  0  =  0.5.  The 
NO-reburning  reaction  significantly  becomes  stronger  as  0 
increases  from  0.5  to  1.5. 

Effects  of  the  chemical  kinetic  scheme 

Although  considerable  investigators  found  that  GRI-Mech  2.11 
is  better  than  GRI-Mech  3.0  in  the  NOx  modeling  66—69],  some 
studies  also  found  that  GRI-Mech  3.0  is  suitable  for  the  MILD 
combustion  modeling  [31,34,38,39,43,75,76].  In  this  section,  we 
examine  the  NO  routes  using  GRI-Mech  3.0  to  investigate  the 
effect  of  the  kinetic  scheme.  The  relative  importance  of  each 
route  is  obtained  by  calculation  with  one  mechanism  at  a  time 
and  the  mechanism  of  each  route  is  shown  in  Table  4.  Fig.  11 
shows  the  calculated  NO  formation  and  reduction  via 
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Fig.  10  —  Calculated  relative  importance  of  NO  formations  and  reductions  via  thermal,  prompt,  N20-intermediate,  NNH,  and 
reburning  routes  of  the  WSR  model  using  GRI-Mech  2.11  mechanism. 
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Fig.  11  -  Calculated  relative  importance  of  NO  formations  and  reductions  via  thermal,  prompt,  N20-intermediate,  NNH,  and 
reburning  routes  of  the  WSR  model  using  GRI-Mech  3.0  mechanism. 


thermal,  prompt,  N20-intermediate,  NNH,  and  reburning 
routes  of  the  WSR  model  of  Cases  5—13  using  GRI-Mech  3.0. 

Comparing  with  ’ig.  10,  it  is  discovered  from  :ig.  11  that  the 
similar  result  is  obtained  by  using  GRI-Mech  3.0,  while  the 
relative  importance  of  the  prompt  route  increases.  There  is 
indeed  difference  of  the  prompt  route  between  two  GRI 
mechanisms.  Lee  et  al.  [66]  found  that  GRI-Mech  2.11  gives 
reasonably  good  agreement  with  the  measured  NO  for  the 
CH4/H2  bluff  body  flame  but  it  under-predicts  NO  levels  in  the 
fuel  rich  region.  They  also  pointed  out  that  GRI-Mech  3.0  gives 
higher  prompt  NO  formation  compared  to  GRI-Mech  2.11. 
Barlow  et  al.  [67]  found  that  GRI-Mech  2.11  provides  decent 
agreement  with  measured  NO  levels  in  lean  and  near- 
stoichiometric  conditions,  but  it  under-predicts  NO  levels  in 
fuel-rich  conditions  of  the  methane  and  air  flame.  They  also 
noticed  that  GRI-Mech  3.0  significantly  over  predicts  the  peak 
NO  levels  for  the  0  =  3.17  and  2.17  flames,  but  it  yields  rela¬ 
tively  good  agreement  with  measurements  in  the  0  =  1.8 
flame.  They  agreed  that  there  is  some  difference  of  the  prompt 
mechanism  between  GRI-Mech  2.11  and  3.0.  Kim  et  al.  [68,69] 
discovered  that  the  GRI-Mech  3.0  significantly  over-predicts 
the  NO  formation  except  in  a  fuel-rich  condition,  while  the 
GRI-Mech  2.11  gives  acceptable  results  on  the  overall.  In  GRI- 
Mech  3.0  the  rate  coefficient  of  reaction  CH  +  N2  N  +  HCN 
is  much  higher  than  that  in  GRI-Mech  2.11.  A  higher  pre¬ 
exponential  factor  for  this  prompt  route  is  the  major  reason 
for  the  over-predicted  NO. 


Conclusions 

The  present  study  has  experimentally  and  numerically 
investigated  the  effect  of  hydrogen  addition  on  the  charac¬ 
teristics  of  MILD  combustion  and  the  NO  mechanism.  When 


the  MILD  combustion  of  NG  and  NG/H2  fuels  is  achieved  in  the 
experiments,  quite  a  uniform  temperature  distribution  and 
very  low  emissions  of  NOx  are  obtained.  The  NO  formation 
and  reduction  mechanisms  under  the  MILD  condition  of  CH4/ 
H2  mixtures  are  examined  with  detailed  chemical  kinetics 
model.  The  main  conclusions  are  summarized  below: 

(1)  In  the  experiments,  with  the  addition  of  hydrogen, 
although  the  furnace  temperature  is  increased,  the 
location  of  the  high  temperature  zone  does  not  change. 
This  high  temperature  region  is  located  at  the  V-shaped 
junction  of  the  upward  fresh  fuel  jet  and  downward 
recirculated  air  jet.  Moreover,  as  the  mass  fraction  of 
hydrogen  in  fuel  (y{j  )  is  increased  from  5.7%  to  14.4%, 
both  the  furnace  average  temperature  and  extracted 
heat  are  increased,  but  the  NOx  emission  remains 
unchanged. 

(2)  From  the  CFD  and  WSR  simulations  it  is  found  that,  for 
lean  conditions  (0  <  0.8),  the  N20-intermediate  route  is 
the  most  important  mechanism  for  producing  NO  under 
MILD  combustion  of  CU.JH2  (taking  60%— 90%  of  the 
total  NO  formation),  regardless  of  the  hydrogen  addition 
(y{j  =  0%— 20%),  initial  mass  fraction  of  oxygen  in 
oxidant  (Yg2  =  3%— 9%),  reactor  temperature 
(Twsr  =  1100  K-1600  K)  and  burner  configuration.  The 
prompt  and  NNH  routes  also  produce  NO  formation  but 
they  are  less  important.  The  NO-reburning  reaction  is 
significant  and  reduces  30%-50%  of  the  total  NO 
formed. 

(3)  For  the  MILD  combustion  at  0  =  0.8,  as  hydrogen  is 
added  to  the  fuel,  the  NO  production  from  NNH  route 
increases  while  that  from  the  prompt  route  decreases. 
However,  the  addition  of  hydrogen  does  not  change 
noticeably  the  NO  production  from  the  N20- 
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intermediate  route.  As  a  result,  as  Y^  is  increased 
from  5.7%  to  14.4%,  the  measured  NO  emission  changes 
little. 

(4)  For  the  MILD  combustion  at  0  =  0.8,  as  the  reactor 
temperature  (TWSR)  is  increased  from  1100  K  to  1600  K, 
the  N20-intermediate  route  remains  dominant  but  its 
importance  reduces  and  also  the  thermal  route  be¬ 
comes  more  important.  The  conversion  from  NO  to  N02 
is  stronger  at  low  temperatures  and  thus  the  relative 
importance  of  N02  to  the  total  NOx  emission  increases 
as  the  reactor  temperature  decreases. 

(5)  For  the  MILD  combustion  at  0  =  0.8,  as  the  initial  mass 
fraction  of  oxygen  is  increased  from  3%  to  9%,  the  N20- 
intermediate  route  increases  its  contribution  in  the  NO 
formation  whereas  those  of  prompt  and  NNH  routes 
decrease. 

(6)  For  the  MILD  combustion  at  0  =  0.5,  the  N20-interme- 
diate  route  is  extremely  important  and  its  contribution 
accounts  for  94%  of  the  total  NO  formation.  However,  as 
0  is  increased,  the  N20-intermediate  route  weakens 
dramatically  while  the  prompt  and  NO-reburning 
routes  become  dominant. 

(7)  The  present  finding  is  insensitive  to  the  kinetic  scheme 
GRI-Mech  2.11  or  GRI-Mech  3.0,  although  the  latter  over¬ 
predicts  the  NO  formation  through  the  prompt  route. 
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